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Abstract : 

By expanding functions of parton fragmentation into a heavy hadron in the inverse of 
the heavy quark mass m,Q we attempt to factorize them into perturbative- and nonperturba- 
tive parts. In our approach the nonperturbative parts can be defined as matrix elements in 
heavy quark effective theory, the shape of the functions is predicted by perturbative QCD. In 
this work we neglect effect at order of ro^ 2 and calculate the perturbative parts at one-loop 
level for heavy quark- and gluon fragmentation. We compare our results from leading log 
approximation with experimental results from e + e~ colliders and find a deviation below or 
at 10% level. Adding effect of higher order in a s it can be expected to reduce the deviation. 
The size of matrix elements appearing at the order we consider for several types of heavy 
hadrons is determined. 



1. Introduction 



Parton fragmentation functions are in general nonperturbative objects in the QCD 
factorization theorem [1] for predictions of inclusive productions of single hadron. For a light 
hadron the fragmentation happens at the energy scale of Aqcd, which is several hundreds 
MeV. Hence it is purely a long-distance process. For the hadron being a quarkonium, which 
consists mainly of a heavy quark Q and its antiquark Q, there is in the fragmentation not 
only long-distance effect but also certain short-distance effect because heavy quarks are 
involved and they provide a large energy-scale, the mass ttlq of heavy quarks. This short- 
distance effect can be well described by perturbative QCD, while the long-distance effect can 
be parameterized with matrix elements of local operators defined in nonrelativistic QCD [2]. 
With these facts an amount of functions for parton fragmentation into a quarkonium is 
calculated (See [3-5] for an incomplete list of references). A question naturally arises that 
can we predict fragmentation function for a hadron containing a single quark Q? Here the 
large mass of the heavy quarks implies certain perturbative effect as it does in the case of a 
quarkonium. In this work we attempt to answer the question. 

Recently our understanding of physics related to hadrons containing single heavy quark 
has grown rapidly. Such achievement is based on the development of an effective theory for 
heavy quarks, i.e. HQET, by starting from QCD (For HQET see reviews in [6]). The basic 
observation leading to HQET is that the heavy quark inside a heavy hadron carries the most 
momentum of the hadron. In this work we will refer heavy hadrons as those containing single 
heavy quark. With this observation one can decompose the momentum of the heavy quark 
into a large component which is roughly the momentum of the heavy hadron, and a small 
component which is at order of Aqcd- By integrating out the dynamical freedom carrying 
the small component HQET is obtained, in which predictions can be expanded in TOq 1 , 
especially, at leading order every heavy hadron has the same mass — tjiq. The difference 
between masses of different heavy hadrons is at order of rriQ 1 . Many applications of HQET 
have been done for weak decays of heavy hadrons. Application of HQET to heavy quark 
fragmentation, i.e., heavy hadron production, appeared first in [7], where an expansion is 
obtained, the expansion parameter is the difference between the masses of the heavy hadron 
and of the heavy quark. The expansion is nonperturbative and formal, it does not tell in 
detail how the fragmentation depends on the energy fraction carried by the heavy hadron. 

For heavy quark fragmentation a heavy quark is the initial parton which is off-shell 
with an invariant mass larger than tuq. One can image the fragmentation as a two-step 
process. Before combining light quarks and glue to form a heavy hadron the heavy quark 
can emit or absorb light quarks and glue. After such emissions it combines light quarks or 
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glue to form the hadron and these light quarks and glue should have momenta at order of 
A-qcd, as the formation is a long-range process. The most momentum of the heavy quark at 
this time will be carried by the hadron and the invariant mass of this heavy quark is at the 
order of the hadron mass. This heavy qaurk can approximately be treated as on-shell. Such 
picture suggests that the fragmentation function may be written in factorized form in which 
a part is for the process from the off-shell- to the on-shell quark and another part is for the 
transition into the heavy hadron. The first part can be treated perturbatively because of 
the large scale which is ijiq at least. The second can be parametrized by using HQET. For 
other initial partons the fragmentation can be thought as a heavy quark is first produced 
and then the formation follows. The production again can be handled perturbatively. 

In this work we use the diagram expansion method to perform the factorization men- 
tioned above. Such method was first used in deeply inelastic scattering to analyse twist-4 
effect [8]. This method can be thought as extension of Wilson's operator expansion to cases 
where the expansion is not applicable, for cases where Wilson's operator expansion is appli- 
cable both methods delivery same results. For readers unfamiliar with this method we refer 
to [8,9,10]. With the method we obtain fragmentation functions as an expansion in TOq 1 , 
in each order the nonperturbative part is contained in matrix elements defined in HQET. 
In this work we will neglect all effect which is suppressed by Wq 2 . For fragmentation func- 
tion D H / a (z), where a stands for the initial parton , H for the heavy hadron and z is the 
momentum fraction carried by H, the expansion may be written as: 

D H/a (z) = D a (z) < 0\O H \0 > +0{\). (1.1) 

m Q 

In Eq.(l.l) Oh is an operator is defined in HQET, its matrix element represents nonper- 
turbative physics. The function D a (z) can be calculated perturbatively. We will calculate 
Dq{z), Dg(z) up to order of a s . The effect at order of tUq 1 is included in the matrix 
element. 

Our work is organized as the following: In Sect. 2 we use the diagram expansion method 
for tree- level diagram and HQET to obtain the factorized form as in Eq.(l.l) for a = Q. In 
doing so, we neglect the difference between masses of heavy hadrons and of heavy quarks. 
It should be be noted that at first look one can keep this difference, but one will have 
some problems at higher order of a s . We will explain the problems in detail. In Sect. 3 we 
proceed to calculate Dq{z) and Dg(z) at order of a s . In Sect. 4 we compare our results with 
experiment at e + e~ colliders and determine the value of several matrix elements < 0|O#|0 > 
by using Z-decays. Sect. 5 is the summary of our work. 

In this work we will use Feynman gauge and <i-dimensional regularization. In this reg- 
ularization infrared(I.R.) singularities appear as poles at d = 4. We take the normalization 
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of a state as: 



<p'\p >=2/(27r) 2 5 3 (p / -p). 



(1.2) 



2. The Factorization and Tree-Level Results for Heavy Quark Fragmentation 

We start in this section our analysis from definitions of fragmentation functions. These 
definitions are first given in [11]. Such definitions are conventionally written in light-cone 
coordinate system. In this coordinate system a d- vector p is expressed as p M = (p + ,p~, Pt) 5 
with p + = (p° + p d ~ l )/\/2, p~ = (p° — p d ~ l )/\/2. We introduce a vector n with n M = 
(0, 1, 0, ■ ■ ■ , 0) = (0,1, Ot) in this system. The function of heavy quark fragmentation into 
a heavy hadron H carries momentum P is defined as[ll]: 

%«W=17 / ^" e "'"" P+ | Tr coio r ^Tr Dirac [n- 7 < 0|Q(0) 

/•OO 

■Pexp{-ig a / dXn-G T (Xn)}a j H (P)a H (P) (2.1) 
Jo 

/•oo 

•Pexp{i# s / dXn- G T (Xn)}Q(x-n)\0>], 

J x~ 

where G^x) = G^(x)T a , G^x) is the gluon field and T a (a = 1, • • • , 8) are the S£/(3)-color 
matrices. The subscript T denotes the transpose. Q(x) stands for the Dirac-field of heavy 
quark. a ! H (P) is the creation operator for the hadron H and _P M = (P + ,P _ ,0t)- For 
hadrons with nonzero spin the summation over the spin is understood. The hadron carries 
a fraction z of the momentum p of the heavy quark as the initial parton, i.e., P + = zp + . 
The definition is a unrenormalized version. Ultraviolet divergences will appear in Djj/q{z) 
and call for renormalization. The renormalization is discussed in [11]. We will use modified 
MS'-scheme. 

At tree-level there is only one diagram in the diagram expansion, which is given in 
Fig.l. We divide this diagram with a horizontal broken line into a upper- and a lower- 
parts. The upper part contains nonperturbative part and we represent it as a black box, 
corresponding to the nonperturbative object Fij(q,P): 

TyfaP) = J d A xe~^ x < 0\Q i (0)a i H (P)a H (P)Q j (x)\0 >, (2.3) 

where i and j stand for Dirac- and color-indices. Because of color-symmetry T(q, P) is 
diagonal in color-space. The contribution of Fig.l can be written as 

Dh/q{z) = ^ J ^^b- nr (Q,P)P^(n-k-n.q). (2.4) 
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If we take if to be a quark Q, the black box becomes a quark line, and F(q, P) becomes 
(27r) 4 5 4 (g — P) (7- q + rriQ), we obtain Dq/q(z) = 5(1 — z). The assumption that the hadron 
carries the most momentum of the heavy quark implies that the dominant ^-dependence of 
the matrix element in Eq.(2.3) is e %x ' p ', the correction to this dependence can be expanded 
in A-qcd/itiq. To proceed further we use HQET to write the field Q(x) as: 

Q(x) = e~ im Q™ ■ {h v (x) + -5-i 7 ■ D T h v (x) + 0(^3-)} (2.5) 

2m Q m z Q 

where 

pn 

D£ = D' t -v Dv», v^ = — (2.6) 

and rriQ is the pole mass of the heavy quark. In Eq.(2.6) is the covariant derivative. 
The effective Lagrangian for the field h v (x) reads: 



It,. „ ,0. „, 1_ 

,2 



L eff = h v iv ■ Dh v + ■^—h v {i'y • D T fh v + 0(-V). (2.7) 



The x-dependence of h v (x) in a matrix element is controlled by the scale at order of 
Aqcd and can be expanded. With Eq.(2.5) and neglecting Mu — rriQ the matrix element 
in Eq.(2.3) can be written as an expansion: 

< Q\Q l {Q)o ] H {P)a H {Q)Q J {x)\Q >= e +%x ' p { < ^\{K)M^ H a H {hv)M\^ > 

- V < 0|(a M / i „),(0)aJ f a H (^) J (0) + (Mi(0)aW(3A>)i(0)|0 > 
1 



+ - < 0|(z 7 • D T h v UQ)cJ H aH{h v )M + {h v )i{<S)<f H a H {il • £> T Mi(0)|0 > 

+ ••■}. 

(2.8) 

The • • • stand for terms which will lead to contributions at order of nig 2 or higher orders. 
The matrix elements in r.h.s. of Eq.(2.8) are matrices in Dirac-indices i and j. The structure 
of these matrices can be determined by using symmetries of parity (P), of time-reversal(T). 
The matrix element in the second and third line can be written as the form: 

< OP^MO^a^yO) + (Mi(O)akatf(0^ w MO)|O > = 

Cg < 0\Tr{((i 1 -D T ) 2 h v )(0)J H a H h v (0)+h.c.}\0>, 
< 0|(z 7 • DThy^a^anihv)^) + (Mi(0)4rM*7 • ArMi(0)|0 > = 

Bij < OlTr^i-f ■ D T h v )(0)a\ja H h v (0) +h.c.}|0 >, 



(2.9) 



where Cg and are matrix elements labeled by Dirac-indices. Because of 7 ■ vh v = h v the 
expectation value in r.h.s. of the second equation is zero. In the first equation we used the 
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equation of motion. The dimension of the expectation value in r.h.s. of the first equation is 
3 in rriQ, while the expectation value in the first line of Eq.(2.8) is 1. Hence the contribution 
from the second line in Eq.(2.8) in final results will be at order of nig 2 . We define the 
operator Oh as: 

Oh = I ^TV{M0)aL(PW(P)M0)}. (2-10) 
With this operator the nonperturbative object T(q, P) becomes: 

F(q, P) = (2n) 4 5 4 (q -P){ 1 -q + m Q )< 0|O H |0 > +0(4^. (2.11) 
Finally we obtain 



m Q 



D h /q{z) = 5(1 -z)< 0|O H |0 > +0{\). (2.12) 

m Q 

In Eq.(2.12) we have neglected all terms which are at order higher than tuq 1 . The obtained 
function is singular at z = 1. If we go to higher orders it becomes more singular, where 
derivatives of 5-function appear. This means that the detail of the shape around z = 1 can 
not be predicted in our approach, but physical predictions still can be made by noting that 
they are convolutions of fragmentation functions with other functions in z, the fragmentation 
functions are distributions. With the fragmentation function given in Eq.(2.12) the integral 
with a test function f(z) is approximated by: 



Jo 



dzf(z)D H/Q (z) = /(l) < 0\O H \0 > +0(\). (2.13) 

m Q 



The situation here may look like the case of inclusive decays of B mesons[12], where one also 
encounters similar expansion as in Eq.(2.12) in which terms from higher orders are more 
singular. In [12] the accurate shape of decay-spectra was interested and the study there is 
corresponding to study the integral in our case: 

I dzf(z)D H/Q (z) (2.14) 

Jl-ZA 

for at order of Aqcd/^q- Because the order of z^ our results in this work can not be 
applied to observables which are related to a integral like that in Eq.(2.14). The observables 
studied later are related to the integral in Eq.(2.13). 

One can perform similar analysis as above for other partons. Based on such analyses we 
propose to write the functions of parton fragmentation into a heavy hadron as the factorized 
form: 

D H/a (z) = t) a (z) < 0\O H \0 > +0(^5-). (2.15) 

m Q 
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In this form the hadron dependence is only contained in the matrix elements, the z- 
dependence is predicted completely by perturbative QCD. It should be noted that the 
matrix element < 0|O#|0 > also contains effect at order of trig 1 , because it is defined in 
HQET in Eq.(2.7) with the accuracy at order of rriQ 1 . In our approach here the function 
D a (z) is just the fragmentation function Dq/ a {z) of a parton a into a heavy quark Q. This 
fact also implies any subtraction for calculating D a (z) at higher order of a s is not needed. 
Usually certain subtractions are needed to extract perturbative parts in a factorized form. 
To examine this we have calculated the matrix element < 0|O#|0 > by taking H = Q in 
HQET at one-loop level. Indeed, the matrix element does not receive any one-loop correc- 
tion. This also shows that the //-dependence of the matrix element is suppressed at least by 
rriQ 2 or by a 2 s . 

In the above approach we have neglect the difference Mh — rriQ in the kinematic as 
in the case with a quarkonium, where the binding energy is neglected in the kinematic. It 
seems that such difference can be kept in the approach. With the difference HQET given 
in Eq.(2.6) is not suitable for our purpose. The reason is: The heavy quark in Fig.l carries 
the momentum q = rriQV + k\. In HQET employed above we neglect k\ at leading order of 
nig 1 . The mass Mh is usually larger than mg, so the fragmentation function is only non 
zero at z = Mh/tuq which is larger than 1. This is in conflict with the definition given in 
Eq.(2.1), which says because of the conservation of momenta that D h /q(z) = for z > 1. 
This problem may be solved at first loo by employing HQET with a residual mass. As 
already noted in [13] that the decomposition of quark momentum q = rriQV + k± where k± 
is the small component, is arbitrary, one can also decompose q as: 

q = (m' Q + e m )v + k[ (2.16) 

where e m is the residual mass at order Aqcd and k[ is the small component. With this 
decomposition one can obtain in analogy to Eq.(2.5): 



Q{X) = e -i(m' Q +e m )v.x . + O(-L)}. (2.17) 

The effective Lagrangian for the field h' v reads: 



L eS = h' v {ivD + t m )h' v + 0{\) (2.18) 

m Q 



Repeating the above steps one gets: 



D H/Q = z 3 5(z-z )<0\O H \0>+O(-^), z = M f . (2.19) 



m Q ' m Q+ e 



rn 
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The operator 0' H is obtained by replacing h v in Oh with h' v . Choosing e m > M~h — m'q the 
function is only nonzero at z = zq < 1. However, the mass m'q is not the pole mass, and 
the pole mass ttiq is the sum m!q + e m . Hence, the choice of e m is not possible. It seems 
that the effect from the difference M H — rriQ can not be handled in perturbative theory. To 
study this effect one may only employ nonperturbative methods or try to sum contributions 
of a series of higher-dimensional operators in [7] . 

In the case of heavy quark distribution in a heavy hadron the difference can be kept 
in our approach. The analysis is similar as that leading to Eq.(2.11), the corresponding 
diagram is just this by reversing Fig.l and p + = zP + . One obtains 

f Q/H (z) = 8{z-^) + 0{X) (2.18) 

where the matrix element corresponding to < 0|O#|0 > equals one plus corrections at order 
of ttiq 2 . The quark line here represents on-shell quark, therefore the problems mentioned 
above will not appear. 



3. Results For Fragmentation at One-Loop Level 



In this section we present a calculation of fragmentation function for a heavy quark into 
a heavy quark and for a gluon into a heavy quark at order of a s . This is the perturbative part 
in the fragmentation into a heavy hadron in Eq.(2.13). In the case of quark fragmentation 
contribution from every diagram at one-loop level contains a I.R. singularity, there is a 
delicate cancellation of the singularity between contributions from different diagrams. We 
show here in detail how this works. 

At one-loop level there are four diagrams contributing to heavy quark fragmentation. 
They are given in Fig.2A-2D. Contribution from each diagram is not only ultraviolet diver- 
gent but also I.R. divergent. However, final result is free from I.R. singularity. With the 
Feynman rule given in [11], the contribution from Fig.2A is: 

Da(z) =—V s / (^) d Tr{ - ig s T a r— ( r, • n(ig s T a n„) 

247r J 27r 7 • (q - k) - m Q + z0+ ^ ^ 

.( 7 . ? + m Q )}.M(n.( P - ? )).^.^-^ 

where fi is the renormalization scale, e = 4 — d, q + = zp + and q 2 = rriq. The contribution 
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from Fig.2B is: 

Db(z) / (~) d Tr{ - ig s T a r — — 7 • n(-ig s T a n„) 

24tt J 2n <y-(q + k)-m Q + z0+ ^ ^ 

■(l-q + m Q )}. - {p _~ Z } _ . Q+ ■ (-1) • 2n5(k 2 )2nS(n • (p - q - k)). 

Both terms are I.R. divergent. In Da(z) the divergence appears at k + ~ 0, while in 
D B {z) the divergence is because the on-shell gluon can carry very small energy. These 
divergences can be regularized in dimensional regularization and are represented by the 
terms as ej 1 , where £/ = d — 4. Performing the loop integration in Eq.(3.1) and Eq.(3.2) 
we obtain: 

= _ |,2, (1 _,).^. r( £)(£). 

■{^-1 + £;+0(e 2 )}, 

£ ' ^ (3.3) 

The T-function with e = 4 — d represents U.V. divergence. The +-prescription is 

as usual. From Eq.(3.3) the sum Da(z) + Db(z) is free from the I.R. pole ej . After 
renormalization the sum is: 

(D A (z)+D B (z)) {R) = l-a s (p){[S(l-z) + -^—} ln£^ + 25(l-,)-2,(M^) + }. 

(3.4) 

The contribution from Fig.2C is just the one-loop correction to Fig.l. Because the 
quark-line is for a on-shell quark, the correction is to external line. This contribution 
contains also I.R. singularity. After renormalization the contribution is: 

D c {z)W = ^5(1 - *) . { A + 7 - H An) - 2 + I In (3.5) 

37T £l 2 fl z 

The last contribution is from Fig. 2D, it reads: 

D d (z) = —/j, £ [ (—) d Tr{(-tg s T a ^) l - 7 • n 7 ^ 

DK ' 24tt^ J K 2n' U ys 1 ^-(q + ty-mQ 1 1 -{q + k)-m Q (3.6) 

• (zs s T^)( 7 • q + m Q )} ■ (-l)(2n)5(k 2 )(2n)5(n .(p-q-k). 
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Performing the /c-integration and renormalization we obtain: 

D D (*)<*) =^M(1 - Z ){ In 4- - 21n(l - z) - 1} 

Oil TUq 

_?^ {[ 2_ Mto)+7] , (1 _, ) + _ii_ } . 

The total contribution to the one-loop correction of -Dq/q is the sum: 2(£)a(z) + + 
-Del- 2 )) +-D_d(-2) ( -' R ' > - In this sum the I.R. divergence in Eq.(3.7) cancells that in Eq.(3.5). 
Therefore the sum is free from I.R. singularity. With these results the function Dq(z) in 
Eq.(2.13) is: 

D Q (z) = D Q/Q (z) = 8(l-z) 

+ -^{[1 - * + - *> + (T^l 4 + " «) (3.8) 
+ 2(1 + ,)ln(l -,) + (!+,) -4(^1^^-^-^-}. 

Now we turn to gluon fragmentation into a heavy hadron. The definition of gluon 
fragmentation function Dh/q{z) can also be found in [11]. Upto the order of ttIq 1 we 
consider there is only one diagram drawn in Fig. 3. It should be pointed out that there are 
more diagrams at higher orders, in which the lower part is connected with the black box not 
only with the quark lines as in Fig. 3 but also with some gluon lines. This is also the case 
for heavy quark fragmentation if we go beyond the order of rriQ 1 . Repeating the procedure 
in Sect. 2 we can obtain 

D H/G {z) = D G (z) < 0\O H \0 > +0{\) (3.9) 

m Q 

as proposed in Eq.(2.13). The function D G (z) is just the fragmentation function Dq/ g {z) 
for a gluon into a heavy quark. From Fig. 3 the contribution reads 

v d-3 



. 2n5(n .(p-q- k)) (-^L^) 2 Tr [ 7 "( 7 • k - m Q )^ ( 7 • q + m Q )] (3 - 10) 

• (p ■ n 9p. P ~ n^(k + q) p )(p ■ ng va - n u {k + q)a)g p ° '■ 
We obtain after integration of the loop momentum k and renormalization: 

D G (z) = D Q/G (z) = -2z + 2z 2 ) In 4". (3.11) 

47T TUq 
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It should be pointed out that the function of heavy quark fragmentation into a heavy 
quark has also been calculated in [14] . In [14] the fragmentation function has been extracted 
by calculating cross-sections at e + e~ collider and the result is formulated in a compact form 
with the -I— prescription. With the property of the +-prescription one can show that the 
result in [14] agrees with that in Eq.(3.18). 

For the later purpose we give here some moments for heavy quark fragmentation into 
a heavy quark. The moment is defined as 

M$\»)= f dzz N - 1 D Q/Q (z^). (3.12) 
Jo 

The moments for heavy quark fragmentation into a heavy hadron can be obtained as 
M^g(/i) =< 0|O#||0 > •Mq V " ) ( / u) in our approach. The first two moments reads 

M^W)=l + 0(a'M), 

m-Vh 1 + ^(>|L-^) + (^)> 

These results will be used in the next section. 



(3.13) 



4. Comparison with experiment at e + e Collider 

In this section we will compare our results with experimental results obtained from e + e~ 
collider. Before confronting to experimental results, we would like to make two comments: 

1) . If the method for factorization used here works for heavy quark fragmentation, it can 
also be applied directly to single heavy hadron production without concept of fragmentation. 
That means that one can expand the inclusive cross-section for production in term of Wq 1 , 
where the same T(q, P) in Eq.(2.3) appears at leading order. In this work we are unable to 
carry out such program. We will still use the results from QCD factorization theorem and 
neglect higher twist effect. In some cases the theoretical analysis will be much simple if one 
uses the results from QCD factorization theorem and fragmentation functions. 

2) . It is confused in the literature about formulation in terms of moments of heavy quark 
fragmentation function for the statement that heavy hadron carries the most momentum of 
the heavy quark. One encounters such formulation for the statement 

Mf = f dzzD H/Q (z) = 1 - O(^). (4.1) 
Jo m Q 
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It is easy to see that the formulation is wrong. For the second moment one can exactly show 
that [11]: 

£<% = !• (4-2) 

H 

With the formulation in Eq.(4.1) the sum-rule in Eq.(4.2) can not be hold. Therefore the 
formulation is wrong. The correct formulation for the statement is 

q^ = l-0(^). (4.3) 

In our approach the correction terms in Eq.(4.3) begin at order of ttiq 2 , and the factor 1 
receives also radiative corrections staring at order of a s . 

If functions for parton fragmentation into a heavy hadron is known or extracted from 
experiment, one can calculate the moments of the functions which is usually at large energy- 
scale. On the other hand one can calculate the moments of the functions using theoretical 
results as those in last sections, where one should take the energy-scale n to be ttiq to avoid 
large logarithmic contribution, and then use the evolution equation to predict the moments 
at the large energy-scale for comparison with the experimental results. Unfortunately, un- 
like parton distributions, the functions are not known well experimentally. However, for a 
comparison one can calculate directly with the theoretical predictions of the functions some 
physical observables, which are well measured in experiment. For this purpose we consider 
the inclusive process 

e+ +e~ -> H + X. (4.4) 

Denoting the beam energy as -Ebeam the variable xh referring to the hadron H is defined 
as: 

x H = = (4.5) 

^beam V * 

where Eh is the energy carried by the hadron H. With the QCD factorization theorem the 
differential cross section can be written: 

da(e + + e~ — > H + X) ^ f dz x H . ^ . , , t . 

— — ' = E J xh T"" { - » )D »>° iz ' "» (45) 

where a stands for all possible partons. In Eq.(4.5) the function h a is perturbative part 
and is known upto one-loop level in MS'-scheme[15]. We will work at one- loop level and 
hence we take only the contribution from the parton fragmentation where the parton a is 
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the heavy quark Q. For the fragmentation function we take the results obtained in the last 
section. The function Kq is[15] : 

M».f) =a«(.){*(l - V) + + -V) ] ^ 

where (Jq(s) is the total cross-section for e + e~ — > Q ± Q at leading order of coupling 
constants in the standard model. The expectation value of any observable 0(xh) as a 
function of xh can be now calculated as 

, Jo dxO(x) J 1 ^j-D h /q(z, li)}iq(xh/ z, fi) 
< 0(x H ) >= 2 1 — 7 

Jo dx fx H -f D H/Q(z, fi)h Q (x H /z, fj) ^ ^ 

Jo dz Jo d y°( z y) D H/Q(z, fj)h Q (y, (j) 
Jo dzD H/Q (z, ft) ■ dyh Q (y, fj) 

If one takes the heavy quark fragmentaion function in last sections and neglect the effect at 
order of Wq 2 , an interesting consequence is that the measured value of observable 0{xh) 
does not depend on the type of heavy hadrons. It should be noted that in experiment one 
can also measure < 0(xq) > by averaging the mean < 0(xh) > for various hadrons H. 
The mean < 0{xh) > is predicted by replacing H with Q in Eq.(4.7). In our approach we 
have: 

< 0(x Q ) >=< 0(x H ) > +0(^5-) =< 0{x H .) > +0(^5") = • • • (4.8) 

m Q m Q 
In experiment, the well studied observable is 0(xh) = %h- Recent measurement at 
^ = M z by ALEPH[16,17] gives: 

< x b >= 0.715 ±0.020, < x Hb >= 0.696 ±0.016 (4.9) 

where H b is the observed hadron in the process (4.4) and it can only be B° or B + . These 
results give certain support for Eq.(4.8). A re-analysis of ARGUS data at ^Js = 10.6GeV 
in [18] also shows: 

< x D >^< x D * >^< x Ac > (4-10) 

at y/s = 10.6GeV. It is interesting to check whether Eq.(4.8) holds for othe type of observ- 
ables in experiment or not. In the following we will concentrate on < xq >. There are two 
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ways to predict < xq >. One way is only to take the perturbative results in Eq.(4.6) and 
in Eq.(3.13) to calculate < xq > in Eq.(4.7). We obtain 



_ m q (v)Jq dyyh Q (y,n) 



< Xq > = < 

m q'W Jo a V n Q{y^) ( 411 j 
a s (u) ( 16 , s 88 1 ^. 

In Eq.(4.9) there is a large logarithmic contribution. This and those at higher orders can 
make the perturbative series unreliable if one only take first two or three orders to make 
numerical predictions. Such logarithmic contributions can however be summed with renor- 
malization group equations. In our case we can take \i = in Eq.(4.7) so that the large 
logarithmic In -j^ in }iq disappears. For Mq(s) in Eq.(4.11) we first calculate them with 
the result in the last section at energy-scale \i = uiq and then use renormalization group 
equation to obtain Mq(s). Here one comment is in order. Since we have already one- 
loop results, one can use renormalization group equations for the moments at two-loop level 
to sum not only the leading log contributions but also next-to-leading log contributions. 
But the equations at two-loop level are unknown. The corresponding equations for par- 
ton distributions are known at two-loop level. At one-loop level there is a simple relation 
between these two sets of equations. At two-loop level this relation is not proven to be 
hold. Therefore we take only the renormalization group equation at one-loop to sum the 
leading log contributions and tree-level results for MQ*\mQ), i.e., we take the leading log 
approximation. The renormalization group equation for Mq (s) reads: 

In Eq.(4.10) 7qq can be expanded in a s (y) and only the leading term is known. The leading 
term for 7qq and 7qq can also be read from Eq.(3.13). We take only the leading term. The 
• • • stands for the contribution from the moments of gluon fragmentation function because 
of operator mixing. Including it the next-to-leading log contributions will be summed. In 
our approach at one-loop level this term should be neglected for consistence. With these in 
mind we obtain 



< x Q >=< x H >= I (l + 0(a,)) (4.13) 

where we used one-loop /3-function for a s and 5 as flavor number. The terms neglected 
at order of a s are those terms: a). The terms at order a s in JiQ(y, ^ = \/s), b). The 
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terms in Mq(h = ttiq) at order a s and c). The next-to- leading term in 7qq , two-loop 
effect in /^-function and effect from the moments of the gluon fragmentation function in 
Eq.(4.12). To make numerical predictions from Eq.(4.13) and Eq.(4.11) we use two-loop 
/3-function for determining a s (ix) at different scales. For li > nib we take 5 as flavor number 
and A^ 5 ) = 200MeV. With these numbers we obtain a s (M z ) = 0.116 which is close to the 
experimental value measured at li = Mz- For li < nib we take 4 as flavor number and 
= 400MeV. With these we obtain a s {m T ) = 0.368 which is also close the experimental 
value at li = m T where m T is the mass of r-lepton. Our input parameters for pole masses 
of heavy quark and for running a s are: 

m h = 5.0GeV, m c = 1.6GeV, A (5) = 200MeV, A (4) = 400MeV. (4.14) 

In experiment there are also data for < xq > measured at y/s = 29GeV. Several groups 
have measured < Xb > and < x c >, where b- and c- quark were identified with their inclusive 
lepton-decays or with charged multiplicity measurements. The results from different groups 
and from different methods are summarized in [19]. We average these results from different 
groups and from different methods by meaning of unconstrained averaging as described in 
[20]. It should be pointed out that these groups have not only measured < xq > but also 
tried to reconstruct the variable z and obtained < zq >. However, such reconstruction relied 
of Monte-Carlo models for fragmentation. We will only make comparison with < xq >. The 
experimental values which we will compare with our predictions are: 

< x b >= 0.715 ± 0.020, <x c > = 0.508 ± 0.011, for Vs = M z , 

< x b >= 0.754 ± 0.034, < x c > = 0.585 ± 0.036, for yfs = 29GeV, (4.15) 

<x c > = 0.640 ± 0.009, for = 10.6GeV, 

where the value for < x c > at y/s = 10.6Gev is from the re-analysis of ARGUS data in 
[18]. The value for < x c > at \fs = Mz is obtained from [21,22] where actually the values 
for < xd* > are measured, we average them and take this value as < x c > according to 
Eq.(4.8). 

With the input parameters in Eq.(4.14) we obtain from the perturbative result in 
Eq.(4.11) the following numbers: 

< x b >= 0.867, < x c > = 0.791, for y/s = M z , 

< x b >= 0.933, < x c > = 0.843, for = 29GeV, (4.16) 

< x c > = 0.906, for = 10.6GeV, 

where fi in Eq.(4.11) was taken to be y/s. Comparing with experimental values there are 
large deviations. The reason is probably because large corrections from higher orders in a s 
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in which large logarithmic contributions exist. Therefore, the predictions above may not be 
reliable. With Eq.(4.13) where the leading log contributions are summed we obtain: 

< x b >= 0.772, < x c > = 0.563, for y/s = M z , 

< x b >= 0.836, <x c > = 0.615, for V* = 29GeV, (4.17) 

<x c > = 0.674, for y/s = 10.6GeV, 

Comparing experimental values for < x b > the deviation from the values given above is 
8% at = M z and 10% at = 29GeV. For the case with c-quark the deviation is 
10% at y^s = M z and is about 5% at other energy scales. Our predictions here are fairly 
in good agreement with experiment. Our predicted values are all larger than experimental 
values. The sources for the deviations mentioned above can be various. Higher twist effect 
neglected in Eq.(4.5) and in Eq.(4.7) can be one of them. An important source is the higher 
order correction in a s . In Eq.(4.13) we neglected this correction. It should be noted that 
the correction in Eq.(4.13) has the form aa s (y / i) + ba s (mo). Because a s {mo) is rather 
large, especially for c-quark, this correction can be large. If we add the corrections from a), 
and b). discussed after Eq.(4.13), the deviation at y^ = M z is reduced to 5% for c-quark 
and to 4% for b-quark. At other energy scales the reduction is not so significant as that 
at y^ = M z , because this correction becomes smaller as y/s decreases. Another possible 
source is the running a s at different energy scale, especially, the value of a s at lower energy 
scales, and also possible nonperturbative effect appearing at these scales for running a s {ji). 
However, a detail study is needed here. 

The last question we will study here is how large is the matrix element < 0|O#|0 > 
defined in Eq.(2.9) for a given hadron. This matrix element should be calculated with 
nonperturbative methods, e.g., with lattice QCD. It can also be extracted from experimental 
results. It should be noted that this matrix element is universal, i.e., it does not depend on 
a specific process. We will use experimental data obtained in Z-decays to extract it for D 
mesons. However, information from experiment is not enough for estimating these matrix 
elements uniquely, certain assumptions must be made. For the estimation we do not use 
the concept of fragmentation. For the inclusive decay 

Z^H C + X (4.18) 

where H c stands for D°, D + and D* + . One can write the branching ratio as 

Br(Z - H c + X) = ^ • P(c - H c ) + ^ • P(b - H c ) + R in (4.19) 

J- z l z 

where we neglected the process of the gluon splitting into cc. The term R{ n is the contribu- 
tion from excited states of H c which are first produced and then decay into H c inclusively. 
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With the method here for the factorization the probability P(c — > H c ) is just the matrix 
element: 

P(c^ H c )=<0\O Hc \0> . (4.20) 

The probability P(b — > H c ) can be expected to be the same as P(c — > H c ) if we consider 
that the b-quark decays first through weak interaction into a c-quark and then the c-quark 
is transmitted into the hadron H c . For D* + the ratio of these two probabilities is extracted 
in experiment which is close to 1[21]: 

P(c^D*) 

— r, ^4 = 1-03 ±0.21. 4.21 

P(b^D*) v ' 

We assume that this ratio is one for all H c . 

For the excited state D* we neglect the contribution from i? in , so we have: 

Br(Z -> D* + X) = Tb ~ b ■ < 0\O D * |0 > . (4.22) 

There is no information for Br(Z — > D*° + X). The matrix element < 0|0_d*o|0 > can not 
be determined with Eq.(4.22). If we assume isospin symmetry for light quarks, the matrix 
element is same as < 0\O D *+\0 >. With the experimental value for Br(Z — > D* + + X) in 
[20] and isospin symmetry we have: 

< 0|O D *o|0 >=< 0|O D «+|0 >^ 0.22. (4.23) 

For D + we take only the decay D* + — > D + + X into account for i? in , according to [20] 
the decay has a chance of 31.7%. With that we have 

Br(Z -> D+ + X) = Tc ~ c ~^ Tbh ■ { < 0\O D+ 10 > +31.7% < 0|O D «+ 10 > }. (4.24) 

1 z 

Taking experimental value we obtain 

< 0|O D+ |0 >^ 0.39. (4.25) 

For D° we take only the contribution for R in from the two decays _D*°, D* + — >■ D° + X. 
The branching ratio for these two decays is 100% and 68.3% respectively [20]. With these 
the branching ration Br(Z — > D° + X) can be written: 

Bx(Z ^D° + X) = TcS ^ Tbl ■ { < 0\O d o\0 > + < 0|O D *o|0 > +68.3% < 0|O D *+|0 > }. 

(4.26) 
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It is interesting to note that with isospin symmetry as assumed before the matrix ele- 
ment with D° is the same as that in Eq.(4.25) and this branching ratio is predictable with 
Eq.(4.26). We obtain 

Br(Z -> D° + X) « 20.1% (4.27) 

which is close to the experimental value 20.7 ± 2.0% [20], or using the experimental value we 
obtain 

< 0|O D o|0 >« 0.40 (4.28) 

which is in consistence with the assumed isospin symmetry. 

Naively one would expect the matrix element <0|Oi)|0>tobe3 times of < 0\Od* |0 > 
because D and D* are spin-0 and spin-1 particles respectively. From our estimation above 
this relation is not hold. The reason for this is that the spin counting can not be applied 
here because the matrix element < O|0#JO > is the probability for the inclusive transition 
c — > H c + X, where the unobserved state X can not be vacuum or a given state and can 
have any possible orbital angular momentum. Further, the unobserved state X for D can 
be different than that for D* . 

There is not data available for b-flavored hadrons, so their matrix elements can not 
be estimated as we did for < 0|Oh c |0 >. However the difference between < 0|Oh c |0 > 
and < 0|OhJ0 >, where Hh stands for B or B* mesons, is at order of m" 1 and of m^ 1 . 
Therefore one can take the value of < 0|Oij c |0 > for the corresponding < 0|O#j0 > as an 
approximation. 

5. Summary 

In this work we studied parton fragmentation into a heavy hadron. We factorized 
the process into a perturbative part and a nonperturbative part. The perturbative part is 
just the parton fragmentation function into a heavy quark at the order we consider. The 
nonperturbative part is a matrix element defined in HQET, which is universal. The z- 
dependence of fragmentation functions is predicted purely by perturbative theory. In this 
work we predicted this dependence for heavy quark- and gluon- fragmentation at one-loop 
level in QCD. With these results we calculated the mean value of the ratio between the 
energy carried by a heavy hadron or a heavy quark and the beam energy at e + e~ colliders. 
Comparing experiment we find that there is a deviation at 10% level between our predictions 
with the leading log approximation and experimental values. The sources for this deviation 
can be several, the important source may be higher order effect in a s as discussed in the last 
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section. The value of the matrix element for D and D* are determined with experimental 
data from Z-decays and one of the branching ratios for Z — > H c + X can be predicted in 
our approach. 

It should be pointed out that our procedure for factorization may directly be applied to 
heavy hadron production with the concept of fragmentation. In this work we do not carry 
out such program for specific process and leave it for future work. 

Note added: 

After the work is finished, the author is informed by Prof. O. Biebel of OPAL group 
about recent measurement of the total branching of c — >• D*. The latest preliminary value 
of the measured branching f(c -> D* + X) is 0.221 ± 0.014 ± 0.013[23]. In the approach of 
our work the branching f(c — > D* + X) is just the matrix element < 0\Od* |0 >. The value 
obtained in this work in Eq.(4.23) is close to the value measured by OPAL. 
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Figure Caption 



Fig.l: The Feynman diagram for contribution to heavy quark fragmentation at tree- level. 
The vertical broken line is the Cutkosky cut, the double line presents the line operator in 
Eq.(2.1). 

Fig. 2 A— 2D: The Feynman diagrams for one- loop contributions to heavy quark fragmenta- 
tion into a heavy quark. 

Fig. 3: The Feynman diagram for contribution to gluon fragmentation into a heavy hadron. 
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